Directed percolation effects emerging from superadditivity of quantum networks 
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Entanglement indcued non-additivity of classical communication capacity in networks consisting 
of quantum channels is considered. Communication lattices consisiting of butterfly-type entangle- 
ment breaking channels augmented, with some probability, by identity channels are analyzed. The 
capacity superadditivity in the network is manifested in directed correlated bond percolation which 
we consider in two flavours: simply directed and randomly oriented. The obtained percolation prop- 
erties show that high capacity information transfer sets in much faster in the regime of superadditive 
communication capacity than otherwise possible. As a byproduct, this sheds light on a new type of 
entanglement based quantum capacity percolation phenomenon. 



Introduction. Percolation (see e.g. [l|, l2|) is a natural 
concept that emerges in the description of spreading pro- 
cesses in the presence of medium imperfections. Percola- 
tion effects in quantum networks have recently been the 
subject of increasing interest [3|-l7[ . Most of the attention 
has been restricted to the generation of large scale net- 
works with maximally entangled states between elemen- 
tary nodes to allow for quantum communication applica- 
tions, starting from initial imperfect, i.e. non-maximally 
entangled state networks. The interesting central new 
insight introduced in [3|, ^ is that local quantum opera- 
tions may be used not only to purify entanglement but, 
simultaneously, to change the topology of the lattice to 
a new one with lower percolation probability threshold. 
This idea has been developed for different states [5|, |6| 
and lattice dimensions [7|,|8|. In a different context, per- 
colation concepts also appeared in quantum information 
theory (QIT) in the study of cluster state generation [9| . 

The capacity of a network determines its utility in the 
domain of communication. The development of QIT has 
led to the uncovering of interesting quantum effects on 
channel capacities, e.g. sup eradditivity of quantum (Q- 
type) channel capacity [lOf. This result followed the in- 
tuition developed in the bound entanglement activation 
effect [11| (where two weak resources activate each other 
becoming collectively useful for some task) continued fur- 
ther in Refs. [1^ and [13|. Independently, the first super- 
additivity effect of classical (C-type) capacity in quantum 
multi-access channels has been described [ij]. Both Q- 
type and C-type supcradditivities have been proven even 
for entanglement breaking channels [l5|. 

In this paper, we consider percolation effects in quan- 
tum networks from a channel perspective. In particular, 
wc show that channel superadditivity can be used to en- 
hance percolation of information through networks. The 
schemes are based on network models of classical infor- 
mation transfer through quantum multipartite channels 
(MACs), which can be mapped to certain types of di- 
rected bond percolation problems [16| . We first consider 
a simple layered communication scheme (A) to demon- 



strate the basic idea of percolation assisted by superad- 
ditive capacities and then describe a more complicated 
scheme of multidirectional communication (B). Interest- 
ingly, the percolation problem in case B docs not seem 
to have been studied elsewhere in the literature. 

The basic ingredients of the quantum networks are: 
a passive fixed underlying network built up of elemen- 
tary entanglement breaking MACs, and an active aux- 
iliary incomplete network consisting of randomly gener- 
ated (open) bonds. The term passive means that no bond 
in the network allows a-priori high capacity communica- 
tion (HCC), whereas elements of the active network are 
high capacity channels. Importantly, the active channels 
can serve to activate HCC through the passive channels. 

Model A - Layered network communication. The ele- 
mentary channel of the passive network is chosen in this 
paper as a 2-sender 1-receiver noisy quantum MAC de- 
picted as the wedge shaped channels in Fig. [TJb): the 
slanted line is a two-qudit (d^-dimensional) input per- 
taining to one user; the vertical line is a single qudit input 
of the second user. The single qudit system is modified 
by one of a chosen set of (orthogonal) d^ unitary oper- 
ations fired by the logical value of the two-qudit system 
(see Ref. 115| for details). This single qudit line is further 
modified by a depolarizing channel and is the sole output 
system at the receiver's end. For moderate to large de- 
polarization - the working regime considered here - this 
channel has poor classical capacity C" for any user |15| : 



C <C Co < Cmax = l0g2 d, 



(1) 



where Cmax is the maximal attainable capacity for a sin- 
gle qudit output, corresponding to an identity channel. 
Communication can be improved by adding a high ca- 
pacity active channel, e.g. an ideal channel, along the 
vertical transmission line. Inputting a two-party maxi- 
mally entangled state to the ideal channel and vertical 
line of a wedge channel increases the capacity along the 
slanted line of the wedge channel to a much higher value 
Co (see Eq.([T])) - a manifestation of superadditivity of 
channel capacities [21| (for high depolarization when the 



a). 



r^^ 



M 



n 




" — ^^r — " — " — y.' 
»i A i ,, ,, ,, 



FIG. 1. Model A - layered communication network: a) passive 
network; users are located at nodes of a square lattice, the lat- 
tice is filled with butterfly shaped primitives, information flow 
is directed from layer t — >■ t-|- 1; b) butterfly primitive consists 
of two MACs (solid and dashed wedges) from Ref. [l5| ; Node 
(user) A^i,* can communicate with A^i-i,t+i, A^i,t+i, A*'i+i,t+i 
with maximal capacity C' <C Co; c) active network filled ran- 
domly by triples of ideal channels allowing HCC (no HCC 
from A to B); d) using entanglement changes the geometry 
of the HCC network leading to directed communication paths 
(black arrows) with capacities > Co from A to B. 
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FIG. 2. Model B - Multidirectional communication networks: 
a) passive network, allowing communication in each direc- 
tion; b) structure of butterfly primitive; c) active channels; 
d) emerging network geometry induced by entanglement be- 
tween active and passive channels. 



channel becomes entanglement breaking, Co « (d+ 1)C" 
[iSlp. We define HCC as transmission at a rate > Cq. 

The structure of the wedge channel induces a sense of 
direction of communication. In the network context, we 
first consider the scenario where only forward communi- 
cation is allowed (see FigdJ . The passive channel inputs 
are shared by nearest neighbour pairs of sites in horizon- 
tal layers giving rise to a butterfly shaped fixed network. 
The active channels are only placed on vertical bonds. 
Suppose that these are available with probability p ~ one 
can assume that identity channels are initially available 
at all vertical bonds, yet due to fragility w.r.t. noise ei- 
ther remain useful (ideal) with probability p or become 
unuseful random operations with probability 1 — p effec- 
tively erasing information. In our scheme, we consider 
a setup where a triple of identity channels is the basic 
active channel 22] ~ this allows the possibility of simul- 



taneous HCC between a given sender and his vertically 
placed receiver as well as between the sender's two hori- 
zontal nearest neighbours and that receiver (Fig. [T]). 

The question of establishing long range HCC in such a 
network is a directed percolation problem. One asks, un- 
der what conditions, is it possible for any user to be able 
to perform directed HCC, through intermediate nodes, 
with a user or users located at a distance scaling with 
the length of the network. One may compare two sce- 
narios - (a) entanglement-free or classical, where no en- 
tanglement is allowed in the protocol used at any node 
and (b) entanglement-assisted (EA), which takes full ad- 



vantage of the superadditive effect described earlier. The 
former case corresponds to communication only along 1~ 
dimensional paths and the percolation threshold prob- 
ability is 1, rendering the network useless for HCC for 
finite loss probability of active channels. The EA scheme 
involves changing the geometry of the HCC network and 
is mapped to a correlated directed bond percolation prob- 
lem where with probability p three directed bonds (form- 
ing an arrow shape) are placed on the lattice (FiglTJe)). 
The threshold probability is significantly suppressed due 
to entanglement induced increased connectivity. We per- 
formed a standard Monte Carlo simulation (see [23|) and 
found the percolation threshold to be Pc = 0.5388 with 
accuracy A = 0.0005. We have checked that the stud- 
ied percolation transition lies, as expected, in the Di- 
rected Percolation (DP) universality class, by comput- 
ing a complete set of critical exponents and found them 
to be in agreement with those obtained for the uncor- 
related directed bond percolation problem. The impor- 
tant basic characteristic of directed percolation is that 
the connected clusters of nodes are geometrically highly 
asymmetrical and are restricted to acute cones (with cone 
angle 7r/2 when p = 1) with axes along the vertical 
lines passing through the starting nodes. Near the crit- 
ical point the clusters are very narrow and essentially 
quasi 1-dimensional, as the probability of obtaining a 
connection with a site a large distance away from the 
source and at an angle 9 from the axis, Q{p,9) > for 
\9\ < 69{p) ^ [p — Pcf' (see e.g. Dhar and Barma [17[). 
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FIG. 3. (a) Comparison of the percolation critical lines 
(p+,pl(p+)) for percolation process on randomly oriented 
square lattice (classical scheme) and butterfly network (EA 
scheme). Dots are numerically obtained points. 



The EA scheme beats the classical scheme in that there 
is a wide window in p for HCC, and that the horizon- 
tal extent of connections \9\ at a distance i — > oo away 
changes from to 7r/4 as p increases from pc to 1. 

Model B - Randomly oriented communication. We now 
move to a general scenario within the described frame- 
work of channels (Fig. ^. Now each node can com- 
municate in both directions with its nearest neighbours 
(main bonds) on a square lattice as well as all of its next 
nearest neighbours (diagonal bonds) through the passive 
network. The configuration leading to this possibility is 
shown in Figl5] (a) and (b) . The active network channels 
are placed along the main bonds of the square lattice and 
are tagged for use in a particular direction. 

HCC between distant nodes is now a multidirectional 
bond percolation problem. We consider the case where 
upward, left-to-right oriented identities are present with 
probability p-^- while downward, right-to-lcft oriented 
identities arc placed with probability p_ . Note, that bidi- 
rectional bonds appear as the independent choice of two 
oppositely directed channels on a given bond. 

This setup is interesting already for the case when, say, 
p+ = 0. Using an entanglement free protocol, capacity 
C > Co can already be obtained through the main bonds 
of the square lattice along the up or right directions, in 
contrast to Model A. This is the well known square lattice 
directed bond percolation problem for which the perco- 
lation threshold Pc ~ 0.64 (see e.g. [17|). But, in the EA 
scheme, due to increased and correlated connectivity, the 



FIG. 4. Comparison of critical exponents r, /3 for classical 
and EA schemes of multidirectional oriented percolation. 



percolation threshold is almost halved w.r.t. the classi- 
cal scheme and has been calculated here to be Pc ~ 0.34. 
Both phase transitions lie in the DP universality class. 
At moderate p w 0.6, the probability of a starting node 
belonging to an infinite or system spanning anisotropic 
cluster (in the thermodynamic limit) Foo > 1/2 using 
superadditivity effects while Foo = otherwise. 

The general multidirectional communication problem 
is the most interesting setup. In the classical scenario, 
this reduces to square lattice randomly oriented perco- 
lation. Problems of this category were studied in the 
context of random resistor diode networks (see e.g. [19[) 
mainly using renormalization group calculations. How- 
ever not many results seem to be available in the litera- 
ture on the p± percolation problem (see however [20[ for 
some analytical properties). 

Here, we map out the phase diagram in the {p^,p^) 
plane in Fig|3] using Monte Carlo simulations (sec Ap- 
pendix). Note that the diagram is symmetrical w.r.t. 
the line p+ = p_ due to system symmetry under the in- 
terchange p_|_ -f^ p_ . This result is to be compared with 
the multidirectional correlated bond percolation phase 
diagram for the EA scheme (Fig. EJd)). As can be seen 
in FiglH the EA scheme is drastically better than the 
classical scheme in the whole parameter plane. 

An interesting feature of both schemes is that they al- 
low switching of universality classes of phase transitions 
from the DP class to the Isotropic Percolation (IP) class 
to which standard bond/site percolation belongs. This 
is facilitated by the choice of parameters p±, and as a 
result one may accordingly change the properties of the 



long range clusters. To provide evidence for this, we cal- 
culate two universal critical exponents (Fig. |4]) /3 and the 
Fisher exponent r (see [23[) as one moves along the criti- 
cal lines (Fig. |3]) of the two models. Recall, that /3 deter- 
mines how the percolation probability Fao ~ [p—PcY i^^" 
creases above the percolation threshold, while the Fisher 
exponent determines how the probability of obtaining a 
cluster larger than size n decreases with n at the critical 
point p = Pc- Fn ^ l/n^^"^. Note that the values of these 
exponents in the DP class are /3 « 0.276, r « 2.112 [ll 
and in IP: P = 5/36 « 0.139, r = 187/91 « 2.0549 [J. 
First consider square lattice multidirectional percolation 
in Fig|4](the classical scheme). Both exponents show that 
as one moves along the critical line, the phase transitions 
lie in the DP universality class until pi ~ 0.4, since this 
region inherits the exponents at the DP point pi = 0. 
At pi = l/2(= p!^), which is an isotropic symmetry 
point, one obtains values of the exponents correspond- 
ing to the IP universality class. In between, there is a 
characteristic crossover region between the two types of 
behaviour. In particular, this means that for a choice of 
parameters approaching the isotropic point, one obtains 
different characteristic growth of clusters (determined by 
different critical exponents) than when in the DP region. 
Secondly the cluster geometrical characteristics change 
from highly anisotropic to isotropic. The butterfly net- 
work percolation problem (EA scheme) basically follows 
the same pattern (see Fig. [4]) and can be considered a 
rescaled version of the classical problem, wherein again 
lies the quantum advantage of the EA scheme. For this 
network, the isotropy point is found to be located at 
P+ = pl ~ 0.225. The slightly lower values of /3 in 
the DP regime as compared to the square lattice prob- 
lem do not seem to be significant - and are a result of the 
extreme sensitivity of the calculated exponents on the 
accuracy of critical probabilities. 

Concluding remarks. We have described percolation 
effects in a channel context showing how directed perco- 
lation effects (not considered before) emerge in the con- 
sideration of quantum networks. In the context of per- 
colation theory, to our knowledge, the multidirectional 
correlated bond percolation problem has not been stud- 
ied before. Finally, quite remarkably, our results provide 
a new entanglement percolation effect in the spirit of Ref. 
[31 . Keeping everything else unchanged, consider the Bell 
measurement (BM) channel of Fig. 1 of Ref. jl5| instead 
of the MAC used here, with AC (BC) playing the role 
of diagonal (resp. vertical) bond of the square lattice. 
Then any randomly generated singlet between B and C 
can be switched, via entanglement swapping, to a singlet 
on the diagonal AC. This leads to the same geometry as 
discussed here but now one asks about the possibility of 
building a long range network of singlets. Note that sin- 
glets are directionless. Without using the BM channel, 
we obtain a "classical" scheme related to square lattice 
percolation which has threshold pc = 0.5. Since direc- 



tionless percolation must certainly be at least as good 
as directed percolation, an EA scheme making use of 
the switching mechanism of the BM channel will have 
threshold at most equal to the calculated threshold at 
the isotropic point of Model B p^ < p+ = pi « 0.225. 
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This is a superadditive effect as the identity is unavailabe 
to the user of the slanted input line of the wedge channel. 

For qubit channels, e.g. this can be physically realized 
as a process of ideal transmission with probability p of a 
single photon having three frequency states and two po- 
larization degrees of freedom for each frequency, where the 
dominating noise process is photon loss. 



[23] See supplemental material for methods. 
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Methods 

The critical percolation properties of the two models were studied using direct Monte Carlo simulations in the 
spirit of Dhar and Barma [l|. We studied the change in behaviour of the probability F„, of appearance of clusters of 
size greater than n, with bond probabilities p and looked for characteristic scaling behaviour expected in the critical 
region to identify percolation thresholds plotted. Below the probability threshold, Fn oc exp(— n) for large n while 
Fn — > constant in the supercritical phase [l|. The critical region is characterized by scaling laws with Fn oc n'^~'^ 
described by a power law dependence at the critical value pc, which we localized by sweeping through super- and 
sub-critical probabilities using an interval bisection method. For model A, cluster size distribution data was obtained 
by performing 10'"' realizations (per value of probability p) of cluster growth starting from a single node. Model B 
simulations were performed on a fixed 2 x 10'^ by 2 x 10^ square lattice also with 10^ realizations for each p, where 
cluster connectivity was identified using a breadth-first search algorithm. 

The qualitative values of the critical exponents t, (3 presented in the paper for model B are defined as follows: 

Fn{p) - n-^^-^'> a.tp = pc 

and 

F^{p)^{p~p,f 

A simple method used to obtain these was to directly calculate the slope of the plots of these two functions in log-log 
scale. Since the problem contains two parameters p+,p_, we chose one of them p_|_ do be the independent parameter 
with p_(]3-|_) determined so as to be on the critical line of the model. For the determination of (3, we considered 
clusters of size n greater than 10^ to be "infinite" or system spanning clusters on the finite lattice. 

Alternatively, we also determined /3 from r and an auxiliary exponent 7, which describes the critical behaviour of 
mean cluster size (n) also readily available from the simulation: 

(n) r^ {pc ~ py^ for p-pc ^ 0^ 

The following universal equation, derived from scaling relations, is known to hold for directed (uncorrelatcd) bond 



inc toilowmg universal equation, derived 1 
percolation [l| and isotropic percolation [2[ 



T-2^ 

■7 



We assumed the equation to be true for the entire critical plane and obtained results for f] in agreement with those 
obtained using the first method. Results obtained in the latter manner are those presented in the paper. 
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